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A theory of tunneling through a quantum dot is presented 
which enables us to study combined effects of Coulomb block- 
ade and discrete energy spectrum of the dot. The expression 
of tunneling current is derived from the Keldysh Green's func- 
tion method, and is shown to automatically satisfy the con- 
servation at DC current of both junctions. 
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Over the past decade, the Coulomb blockade of tun- 
neling has attracted growing interest mainly because of 
its possible applications in single-electron electronics jj] . 
In view of ongoing advances in microfabrication tech- 
niques of semiconductor devices whose discrete energy- 
level structures are prominent, it seems urgent to de- 
velop a general theory of tunneling which allows us to 
analyze combined effects of Coulomb blockade and dis- 
crete energy-level spectrum of the quantum dot. While 
a number of articles study the combined effects |^-^| , in 
a usual perturbative treatment, the current conservation 
has been satisfied by adjusting a fitting parameter such 
as a electrochemical potential of the island. By general- 
izing the Keldysh equation developed by Caroli et al. , 
we obtain an expression for the I-V characteristics such 
that the currents at the left and right junctions satisfy 
the conservation of DC current j|] . We also find that the 
half-width of the resonant peak( written by T nn i below ) 
is strongly modulated by the Coulomb interaction. 

The system consists of two tunnel junctions with ca- 
pacitances Cl and Cr, and a gate capacitor C g con- 
nected to the central island. The Hamiltonian consists 
of the electronic part He, the transfer part Ht, and the 
part Hem that describes the electromagnetic (EM) envi- 
ronment surrounding the junctions. The electronic part 
consists of electrode and island parts, 



He — -^kaC kQ CkQ- 



mn 
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where a represents a set of parameters that, together 
with wave vectors k, completely specify the electronic 
state of the left (L) or right (R) electrode, and m, n spec- 
ify the energy levels of the central island. The symbol 
S mn denotes the Kronecker's delta and describe ma- 
trix elements of impurity scatterings in the island. The 
transfer part is described by 



H^ — 



[V n] < a (t)cl a e l ^d n +h.. 



where <j) a (a — L, R) are the phase operators which 
are canonically conjugate to the charge operators Q a 
on the junctions, satisfying the commutation relations 
[4>a, Qp] = "ieSap ||- The EM part is described by 
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+ H({4> a }) +const., (1) 



where C = C l Cr/(C l + Cr), Cx = C L +C R +C g , C x = 
C^Cq /[Cl(Cr+C 3 )], and H({cf) a }) describes the phase 
part of the Hamiltonian whose concrete form depends on 
the external circuit. In the second line of Eq. ([!]), new 
operators q and Q are introduced such that 

Ql = Cr + C g Q g - C g V g 



}r_C l Q q~C g V g 



Cr 



Cs Co 



Qg _ CL Q <1 ~ CgVg 



C, 



Ce Co 



Cy 



v„ 



The commutation relations for the new operators are 
given by [if), q] = ie and [tp, Q] = ie, where 



<t>L = lp + K L tp, 



-if) + KRlfi 



(2) 



with k l = C R /{C L + Cr) and k r = C L /(C L + Cr). 
The current at junction a (=L or R) is given by [fLOf 



J a (t) = {-If | Re|^F„ ka (t)G< ka (t,t)| , 

I kri ) 
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where (3=0 for the left junction and f3=l for the right 

junction and G< kQ (M') = «(cL( t ) e_i *° (t) d n {t')) is an 
analytic continuation of the contour-ordered Green's 
function GjX (r, r') which is defined in the interaction 
representation by 



G nkQ (r,r') ee i(T c {c{ a {T')e~^ T ">d n {T) 



x exp 
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where Tq is the contour-ordering operator. We assume 
that electrons in the left and right electrodes are nonin- 
teracting. Then the only nonvanishing terms in Eq. (^) 
are those in which cl (r') is contracted with c kQ ,(ri) in 
the exponential term. 

GnUr,r') = [ ^E^k Q (r 2 )<Tc{cL(r)c kQ (T 2 )}) 

x (T c {e^Vj^d n (T)dl(r 2 ) exp(-^^H T (n)dT^ }). 

To proceed with calculation, we decouple the term 
(T c {e- 4 ^e i ^ r >>}) on the right hand side of Eq. (g). 
In spite of this decoupling, the conservation of DC cur- 
rent is satisfied as will be below. We then obtain, after 
analytic continuation to real time, 



G< kQ (M') = E 



dtx 



(6) 



where G< m (t,t') = i{d\ n {t')d n {t)) is the Green's 
function of the central island, and <? kQ ,(^i, £2) and 
9k a {h,h) are given in terms of g£ a (ti,t 2 ) = 
i&MZkoih)) and P<{t x M) = (e^^e-^M) 
as 9L(ti,t 2 ) = gL(tiM)P<{tiM) and ~gl a {tx,t 2 ) = 
^2 - hM^M^ihM) - 9ii a {txM)P>{h,t 2 )]- 

Here g^ a (h 7 t 2 ) = -i(cka(h)cl a (h)) and P>{t l7 t 2 ) = 
( e -^(*i) e A,(ij)), Substituting Eq. (§) into Eq. (|), we 
obtain 

/oo 
dt'V* ka (t'){G> m (t,t') 
- OO 

Kmn 

X9L(t',t)P<(t',t) - G< m (t,t')g> a (t',t)P>(t',t)}, (7) 

where G> m (t,t') = -i(d n {t)d\ n {t')) . By introducing 
a noninteracting self-energy, Eq"", defined by Gq < = 
GgE^Gg, and using the Dyson equations, (1 + 
G r E r )Go = G r and G£(l + E a G a ) = G a , we may cast the 
Keldysh equation G>< = (1 + G r IT)G5 ,< (l + E Q G a ) + 
G r E><G a into the following form Q: 

Gnmitjt') = E / ^*l^*2G^ ril (t,tl)E^ nin2 (t 1 ,t 2 ) 

xG° 2m (t 2 ,0, (8) 

where 

(9) 
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The first term on the right-hand side describes the self- 
energy of the island in the absence of disorder, interac- 
tion, and tunneling. Since it corresponds to the free part, 
it is infinitesimal: E>„<„ 2 (e) = 2W m „ 2 r?[/o(e)-l/2=Fl/2], 



where 77 is a positive infinitesimal, > (<) refers to the 
minus (plus) sign, /o(e) = (exp/3(e — Vd) + and 
Vd = (ClV + C g Vg)/Cz. The second term describes a 
self-energy due to tunneling: 



S Tnin 2 (*l>*2) 



E 



v* ka (h)v n2]ca (t 2 ) 



g^(h,t 2 )P><(h,t 2 ), (10) 



(^)and the third term describes effects of scattering: 
Zf< in ,(e) = iSn in . 2 [f d (e) - 1/2 =p 1/2]/t b , where r s is a 
scattering time in the island and /d(e) is a Fermi distribu- 
tion function in the presence of the scattering. Because 
of the infinitesimal factor rj, the free part is important 
only when the remaining parts are absent. 

In the following we focus on the case in which 
r nnlka (t,t') = {2it/n)V^ ka {t)V n >i ta ^ is a real func- 
tion of t - t'. From Eq. @), J L - J R = 
e Ho ^'( G n™(M0^ m « - G< m (t,t')Z> mn (t>,t)). 



Because G><(t,t') contains E>< m (i,t') (Eq. (JT0|) ) , the 
conservation of current through the two junctions is auto- 
matically satisfied, that is, Jl = Jr- This is a great step 
forward, considering the fact that in a usual perturbative 
treatment one has to adjust the electrochemical potential 
of the island so as to satisfy the current conservation. 

The retarded and advanced Fourier-transformed 
Green's functions at the central island, G r nn , (e) and 
G^ n /(e), are derived from the Dyson equation: 



HGZ^T 1 = %£ a (e)r 1 -as; 
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where T nn >(e) = 2ImE| rm , (e), and 2nhT,j, nn , (e) = 
Ska r rm'ka(e)<?kQ(e). The real part of the self-energy 
A n „/(e), which are due to scattering and tunneling, shifts 
energy levels in the central island and we will regard this 
effect as included in our assumed one-body energy lev- 
els of the quantum dot. Here T nn /(e) broadens resonant 
peaks and is strongly modulated by the Coulomb inter- 
action through the P y< functions: 

/d 
^r n „, kQ ( ei )[(i-/ Q (£ kQ )) 

x P>(e - a - £ ka ) + f a {E^ a )P<{e - ei - E ka )], (12) 

where / L (e)=l/(e /3 ( £ -- E ^ y ifl) and f R (e)=l/(e^ e - Ep h-l) 
with Ep being the Fermi energy of the electrodes. 

We show here that the free part of the self-energy can 
be neglected even in the absence of scattering. Because 
Gl^ n = G\Y><G% is given by 



Gl 



(e-E n -V d +A n (e)) 2 + (r, + r n (e)/4y 



(13) 



^free n( e ) remains nonvanishing only when r„(e) = as 
can be seen from the relation, lim.^+o r]/[(e — E„ — Vd + 



2 



A„(e)) 2 + V 2 } = *5(e - E n V d + AJe)). If E><(t) = 
Eka r nka(05ka < (^) p a >< (i)/27r, Eq. (@) is cast into 

^=(-l)^^/^{G>(ei)E< a (e 1 )-G<(e 1 )S> a (e 1 )}. 



(14) 
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Because r„(e) = J2 a C^na( e ) ~ ^ Q (e)) = means that 
£> Q (e) = £< Q (e) = [note that each term in Eq. |H| 
is non- negative] , the current disappears when r„(e) = 0. 
Thus the free part can be neglected in the expression of 
current. In a usual resonant tunneling problem, T„ = 
+ T R (width at half maximum: constant), so rj in the 
denominator can be neglected anyway. 

With A nin2 (e) ee i(G r nin2 (e) - G a nirl2 (e)), and 
he nm (e) = G^ ni (e)G^ 2m (e)/A nin2 (e), J L in the absence 
of scattering can be cast into the following form, 



Jl 



2nh 



E 



KT ki (e 1 )vw lki (e 1 )y n * ik , i? ( ei )y n2k 



kk'nmni ri2 

x {h(E kL ) (f R (Ev R ) - l)P^(ei - E kL )P>(e 1 - E k , R ) 

This is the central result of this paper. In the absence 
of the charging effects, i.e., Pa < i e — €> ) = $( € — Jl 
is reduced to the familiar expression of resonant tunnel- 
ing §J). When S^ Q „(e) = 5 mn E% a , all Green's functions 
are diagonalized. This corresponds, e.g., to a situation 
in which energy levels in the island are mutually uncor- 
related during the tunneling process. Then A nin2 (e) re- 
duces to 



A n (e) 



nr n (e) 



[e-(E n + V d )-A n (e)] 2 + [T n (e)] 2 /4 



(16) 



and 6„(e) reduces to r„(e) 1 . 

From Eq. (§), we see that P> can be divided 
into two parts: P>(t 1 ,t 2 ) = ( e - lK "^ tl '> e iK "^ t2 '>) 
x ( e -< Q ^(ti) e < Q Vi(i2)) (£ £ = i an d £ R = wnere t he 
first factor describes influence of the electrodynamic en- 
vironment (i.e. the external circuit) j^], and the sec- 
ond factor describes quantum fluctuations of the is- 
land charge. Using the relation {m\e~ 11 ^^ e ll ^^\m) = 
e -if [i+2 7 (m-n c )] tj where fi c =(C L V+C g Vg)/e, 7 is an ar- 
bitrary constant and \m) is the charge eigenstate of 
the island with charge em, we obtain 

p -/3U(m-n c ) 2 

( e -i-Y>P(i) gnVKO)) = y^_ [l+2 7 (m-n )]t 



-e 



where U = e 2 /2Cs and c u ee EHe^" 1 " 1 "* 1 ' 2 , and [m] 
runs over 1, 2, ■ • ■ , 2iV, where N is the number of doubly 
degenerate energy levels. Fourier transforming this in the 
high-impedance limit gives 



-0U(m~n c ) 2 



-2nh6(eTE><±K 2 a E C!C ), (17) 



where E>< = U(l ± 2[m - n c }( a ), and E Cm = e 2 /2C x . 

Before we discuss the main consequences of Eq. (|l5|), 
let us discuss the relationship of our theory to the 'ortho- 
dox theory' of double junctions (lj] . In the orthodox the- 
ory, the tunneling rate is assumed to be small, and tun- 
neling processes at the left and right junctions are treated 
separately by assuming that electrons on the island are 
always at thermal equilibrium. To ensure this, we as- 
sume that inelastic scattering in the island is large and 
that tunneling is weak so that 2TT/h\V* ka (e)V*, ka (e)\ < 
%/2t e <C Ep, where Ep is the Fermi energy of the cen- 
tral electrode. Equation (Q) can then be shown to reduce 
to the corresponding formula in the orthodox theory. In 
fact, under these assumptions, we have Im[G^' a (e)] _1 ~ 
±h/2r s and G >< may be approximated as Gq < because 
neare ~ Ep we have lim n/Ts ^ +0 h/ (2t s ) /[(e- E n -V d ) 2 + 
(?i/2t s ) 2 ] = 7r<5(e— E n — Vd). Thus, with appropriate iden- 
tification of parameters, our formulation reduces 

Amn2 (ei)6nm(ei) 
(h(E kL ) - l)f R {E^ R )P>{e 1 -E kL )P<(e 1 - E k , R )} . (15) 

to that of Ref. In the following discussions, we 

consider the case in which the self-energy is diagonal [i. e. 
we assume Eq. (|l6|)] and r ml 'ka(e) does not depend on 
n,n' and k, namely, T nn : ka (e) — r Q (e) and C g — > and 
V d = C L /Cj:V + V g . 

Figures |l| and || show I-V characteristics and the I-V g 
characteristics, when the charging energy U = e 2 /2Gs 
is much greater than the discrete energy-level spacing 
AE n . From the relation T Q (e) = 2TrD k (e)\V ka \ 2 , resis- 
tances R a (a — L,R) can be evaluated to be R u /Rk = 
(T a (Ep)Dd(Ep))~ 1 , where Rk is the resistance quan- 
tum h/e 2 =25.8kft and D d (E F ) is the DOS of the island 
at the Fermi energy. The resistances are estimated to be 
R R /Rk ~ 40(weak tunneling) when D d (Ep)Ep ~ 1. 

Figure [| shows the I-V characteristics in which AE n 
=25meV larger than the elementary charging energy 
C/=0 or lOmeV. Here we assume r Q (e) = ^e/EpTao 
for the energy dependence of the density of state (DOS) 
Dk(e) of the electrodes on the I-V characteristics. In 
the absence of the Coulomb interaction, i.e., U = 0, 
equally-spaced resonant tunneling peaks are seen. The 
Coulomb interaction deforms this I-V characteristic in 
a manner depending on the strength of the interaction 
and temperature. As the Coulomb interaction increases, 
the deformation of the peaks becomes more pronounced, 
and some peaks are suppressed by large Coulomb gaps. 
At T=100K with U—lOmeV, the resonant peaks are 
thermally blurred. The dips at 0.2 eV and 0.4 eV ap- 
pear where original resonant peaks overlap with Coulomb 
gaps. Our results shows that the complex I-V character- 
istics appear due to combined effects of Coulomb block- 
ade and resonant tunneling in the quantum dot. 

In conclusion, we have studied combined effects of 
Coulomb blockade and resonant tunneling through a 
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double-barrier system by using the Keldysh equation. A 
general expression of the I-V characteristics is obtained 
which automatically satisfies the conservation of DC cur- 
rents at both junctions. We found that, in a resonant- 
tunneling-dominated regime, some of resonant tunneling 
peaks are suppressed due to the charging effects. 
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FIG. 2. Coulomb oscillations in the continuum limit for 
r.L=5meV, r.R=0.1meV with varying values of the elementary 
charging energy U = e 2 /2C s , where, E F =0.2eV, T=4.2K, 
and V=0.01eV. 

FIG. 3. Combined effects of the Coulomb blockade and 
resonant tunneling for T = 1 and 100K, where Ef 
=0.07eV (T=1K), V a =-0.05eV, r L0 =6.0 x l(T 4 meV, and 
r H0 =4.0 x lCT 4 meV. The curves at T=100 K refer to the 
right scale and those at 1 K refer to the left one. 
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FIG. 1. I-V characteristics in the continuum limit for two 
different gate voltages V g =0 and -0.5eV and two temperatures 
T=4.2 and 30K, where £ F =0.2eV, U =5meV, r L =5meV, 
r H =0.1 meV. 
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